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Abstract 

In a cooperative multiple-antenna downlink cellular network, maximization of a concave function of user rates 
is considered. A new linear precoding technique called soft interference nulling (SIN) is proposed, which performs 
at least as well as zero-forcing (ZF) beamforming. All base stations share channel state information, but each user's 
I message is only routed to those that participate in the user's coordination cluster SIN precoding is particularly 

' useful when clusters of limited sizes overlap in the network, in which case traditional techniques such as dirty paper 

coding or ZF do not directly apply. The SIN precoder is computed by solving a sequence of convex optimization 
problems. SIN under partial network coordination can outperform ZF under full network coordination at moderate 
SNRs. Under overlapping coordination clusters, SIN precoding achieves considerably higher throughput compared 
' to myopic ZF, especially when the clusters are large. 

Index Terms 

Cellular network, convex optimization, cooperation, interference mitigation, linear precoding, multiple-antenna. 
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I. Introduction 



> ■ Interference management is a fundamental challenge in wireless cellular systems. In this paper, we 
consider the downlink cellular network, and investigate the performance benefits of allowing cooperation 
pr^ and joint processing among the base stations. Without base station cooperation, the system is interference- 
CN limited, i.e., the signal-to-interference-plus-noise ratio (SINR) at the mobiles cannot be improved simply by 
op increasing the base station transmit power, since higher transmit power also creates stronger interference. 
Given the deployment of a fixed number of base stations and antennas, one approach to increase system 
throughput is to allow the joint encoding of user signals across the base stations. In this case, assuming 
perfect cooperation among the base stations, the downlink system can be modeled as a broadcast channel 
(BC). However, the theoretically optimal dirty paper coding (DPC) transmission scheme for the BC can 
be too complex for practical implementation. Zero-forcing (ZF) beamforming is a simple linear precoding 
5^ , technique that offers good performance in a BC. In this paper, we propose a new linear precoding technique 
called soft interference nulling (SIN) that performs better than or equal to ZF. The SIN precoder can be 
found by solving a sequence convex optimization problems. The SIN precoding technique applies to 
the case when the terminals have multiple antennas, as well as the case when the users are served by 
overlapping coordination clusters in the network, where each coordination cluster consists of a limited 
number of cooperating base stations. SIN precoding is particularly useful in the setting of overlapping 
clusters, since under this scenario traditional precoding techniques such as DPC or ZF cannot be applied 
directly. 

In wireless communications, multiple-input multiple-output (MIMO) transmission techniques have been 
shown to provide substantial improvement in channel capacity [[T]|, |l2l. For cellular downlink networks, the 
system throughput using time division multiplexing access (TDMA), ZF, and DPC are compared in [[3]|, 
im, and the performance of ZF is studied in [[5]|, [[61 . Moreover, ZF is generalized for multi-user MIMO 

The material in this paper was presented in part at the Forty-Seventh Annual Allerton Conference on Communication, Control, and 
Computing, in Monticello, XL, September/October, 2009. 
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channels in Q, |[8l. Different preceding schemes for MIMO BCs are presented in [|9l, [fTOll . The optimality 
of DPC in a MIMO BC is shown in [ITTII . For single-cell multiuser MIMO channels, the optimization of 
different performance metrics in terms of the user rates or SINRs are considered in fT2l- |fT6l . In [[T7], the 
transmit precoder is designed to maximize signal strength relative to the interference it causes. Cooperating 
base stations with overlapping coordination clusters in the cellular uplink channel is considered in [fTSl . 
When the user signals are jointly encoded by separate base stations, they are under per-antenna power 
constraints (PAPC). ZF under PAPC are considered in |191-|121j|, and DPC under PAPC is treated in [|2l. 
System-level performance gains in collaborative networks under full-network coordination are studied in 



We consider the maximization of a general concave utility function of user rates in a cellular downlink 
network. We focus on linear precoding techniques, under the assumption that interference is treated 
as noise. We assume all base stations share channel state information (CSI), but a user's message is 
only routed to the base stations that participate in the user's coordination cluster. The required backhaul 
bandwidth to share CSI between the base stations depends on the mobile speeds and update frequencies. 
For typical applications, the requirement for sharing CSI is much less than that for sharing user data. In 
particular, in [25 J, it is shown that the requirement on backhaul bandwidth is about an order of magnitude 
less for pedestrian speeds. 

The remainder of this paper is organized as follows. The MIMO cellular network downlink model 
is presented in Section Ull Section |lll] explains base station cooperation and the precoder optimization 
framework. Section |W] describes the soft interference nulling algorithm and different clustering tech- 
niques. Numerical results on the performance of different downlink precoding algorithms are presented 
in Section |Vl under the settings of a line network and a cellular network. Section |VI] concludes the paper. 

Notation: In this paper, (R+) R is the set of (nonnegative) real numbers, C is the complex field, 
and 1 denotes the two-element set {0, 1}. Dimensions of vectors/matrices are indicated by superscripts. 

(H^^) is the set of N x N positive semidefinite (definite) Hermitian matrices. 1^ is the N x N 
identity matrix; Omxn is an M x N matrix with all zero entries. A^, , are the transpose, conjugate 
transpose, and pseudoinverse, respectively, of a matrix A. [A]ij is the matrix's entry, and [aij] refers 
to the matrix comprising the entries Uij. The matrix diag(a) is diagonal, with its diagonal given by the 
vector a. The operators E[ ■ ], det, tr denote, respectively, expectation, determinant, and trace. For random 
variables, x ~ CJ\f{fi,Q), where G C^, Q E H^, means that x is a circularly symmetric complex 
Gaussian random A^-vector about mean n with covariance matrix Q. 



Consider a MIMO downlink cellular network as depicted in Fig. \T\ Suppose there are B base stations 
and K mobile users in the network. Base j has Mj transmit antennas, j = 1, . . . , B; and User i has Ni 
receiver antennas, i = 1,. . .K. We consider a narrow-band fiat-fading channel model. Wireless systems 
with wider bandwidth may be modeled as multiple narrow-band channels using modulation schemes 
such as orthogonal frequency-division multiplexing (OFDM), and most techniques discussed in this paper 
remain applicable. The MIMO complex baseband channel gain from Base j to User i is denoted by the 
matrix Hij E C^'^*^^. Suppose Xj E C^^^ is the transmit signal at Base j, and yi E C^* is the receive 
signal at User i, then the discrete-time downlink channel is described by 



where each Zi ~ CJ\f{0,lNi) E C^' is independent zero-mean circularly symmetric complex Gaussian 
(ZMCSCG) noise. 

We consider a block-fading channel model: the channel gains realize independently according to their 
distribution at the beginning of each fading block, and they remain unchanged within the duration of 
the fading block. In this paper, we assume the channel states can be estimated accurately and conveyed 



m, f24j. 



II. System Model 
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Fig. 1. Downlink cooperative MIMO cellular network with limited coordination clusters. Base j, with Mj antennas, is under power 
constraint P,; User i, with A'^; antennas, has channel Hij from Base j. In the figure, the coordination clusters are: B\ = B2 = {lj2}, 
Bi = {2, 3}, and the corresponding user sets are: Ki = {1, 2}, IC2 = {1, 2, 3}, /C3 = {3}. 



in a timely manner to all base stations: i.e., the channels are known at all terminals. Each Base j is 
under a transmit power constraint of Pj. We consider a short-term power constraint: i.e., E[a;^Xj] < Pj, 
j = 1, . . . , B, where the expectation is over repeated channel uses within a fading block; power allocation 
across fading blocks is not considered. We assume each fading block is sufficiently long so the transmitters 
may code at near channel capacity using random Gaussian codewords. 

III. Cooperative Cellular Networks 

A. Base Station Coordination Clusters 

Let us consider the scenario where a mobile user is served by a cluster of cooperative base stations. 
Particularly, each User i specifies a coordination cluster Bi, where Bi C {1, . . . , B} is the set of base 
stations that participate in the cooperative transmission to User i, i = 1,...,K. Different clustering 
techniques are discussed subsequently in Section |TV-B[ We assume all base stations share global CSI, i.e., 
all bases have knowledge of all channels Hij. Furthermore, the base stations in the coordination cluster Bi 
all know the message intended for User i, and they may jointly encode the message in their transmission. 
Conversely, the base stations unaffiliated with User i's coordination cluster {j \ j ^ Bi} do not have access 
to User i's message. Let the total number of transmit antennas in coordination cluster Bi be denoted by 

M, ^^M,-, z = l,...,K. (2) 

We use the notation Bi[n] to represent the nth element in the set Bi (sorted ascendingly), and \Bi\ to 
denote the set's cardinality. Note that the coordination clusters for different users may overlap, i.e.. Base j 
may participate in the transmission to multiple users. We denote the set of users served by Base j as 

}Cj^{z\j eB,}C{l,...,K}, j = l,...,B. (3) 

For example, the coordination clusters of Users 1,2,3, and their corresponding user sets, are illustrated 
in Fig. [B 

B. Linear Precoding 

In this paper, we consider linear precoding at each base station. Specifically, let the transmit signal at 
Base j be given as follows: 

Xj = ^ GjiUi (4) 
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where Ui E C^^' represents the information signal from coordination cluster Bi to User i, and Gji E C*^^ 
is the precoding matrix for User i's signal at Base j. Note that in the above formulation, we allow each 
User i's information signal to have multiple spatial streams (up to the total number of transmit antennas 
Mi in its coordination cluster). Without loss of generality, entries of a user's precoding matrix may be 
set to zero if the spatial streams are not all active. Such representation is particularly useful under limited 
coordination clusters: since the clusters may overlap, it may not be clear a priori how many spatial 
streams are supported for each user, and how many antennas are used for beamforming vs. nulling out 
interference. We assume each spatial stream consists of independent data signals, i.e., E[uiuf^] = If^^, 
and E[uiU^] = Oaj-xm^ for i 7^ fc. It is convenient to specify the design of the precoding matrices Gji in 
terms of their corresponding covariance matrices 



Qi ^ G,Gf E H 



G, 



G 



E C 



MiXMi 



(5) 



where i = 1,. . .,K, and Bi[m] is the index of the mth base station in coordination cluster Bi, with 
m = 1, . . . ,\Bi\. To generate transmit signals with the specified covariances, we may set the precoding 
matrices G, to be 



G, = K A 



-/2 



K 



(6) 



where Di is diagonal, and ViDiV/^ = Qi is the eigendecomposition of Qi. 

For each User i, the desired signal is J2jeB HijGjiUi, and the interference from the other users' signals 
is given by J2k=i kjti J2jeBk ^ij^jk^k- Note in the above that each user suffers interference from all base 
stations: i.e., a frequency reuse factor of 1 is assumed. In this paper, we consider the design of linear 
precoders where interference is treated as noise (i.e., interference pre-subtraction schemes such as dirty 
paper coding are not considered). Let Ri be the achievable rate for User i. When other users' signals are 
treated as noise, the following rate is achievable |l2l, |[26ll using Gaussian signals 



det(/jv, + Ef =i H.GkQkGj^H, 



Ri = log ■ 



det{lN. + j:t=i,k^.H,G,Q,G^Hn 
where Hi is the aggregate channel matrix from all base stations to User i 



K 



(7) 



Hi = \Hii Hi2 



HiB] E C 



N,xM 



(8) 



and Ck are constant matrices (each entry of is either a or 1) that represent the association between 
the base stations and the users in the coordination clusters 



Gi = [Eb,ii] 



l3^[2] 



where 
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The M X M matrix CiQiCf represents the full covariance matrix of User i's signal with respect to the 
transmit antennas of all base stations in the network. Effectively, the matrix Ci stipulates that the non- 
participating base stations have precoding weights of zero for User i's signal. In terms of the matrices Ci 
and Ej, the transmit power constraint for Base j can be written as 



K 



Y,^v{EfCiQ,CjE, 



< P. 



3 



B. 



(11) 



To illustrate the construction of the association matrices Ci, consider, for example, a cellular network with 
B = 5 base stations, with each base station having a single antenna. Suppose the coordination cluster for 
User i is Bi = {1,3,4}. Then User z's signal is characterized by the covariance matrix Qi E H'^, which 
describes the joint signal from base stations 1,3,4. The association matrix Ci and the full covariance 
matrix CiQiCf, respectively, are 



CiQiC- 



■T 





[Qiki 













mi.2 



m2,2 





mi,3 



m2,3 





(12) 



Given the aggregate channel matrices, association matrices, and covariance matrices {Hk, Ck, Qk, for 
k = 1, . . . , K), the MIMO rate Ri in (|7]) can be achieved through singular value decomposition (SVD) 
BH, or minimum mean square error (MMSE) detection with successive interference cancellation (SIC) 
[[ll. The SVD and MMSE-SIC strategies are capacity-achieving (when interference is treated as noise) 
for arbitrary numbers of transmit/receiver antennas, and number of spatial streams being active. 



C. Optimal Precoder 

Let the user rates in © be denoted by the vector R = \Ri . . . RkY ^ We are interested in 
maximizing a concave utility function of the user rates. Let U : R denote the utility function, 

where U is concave on R^. Concave utility functions can be used to model a wide class of resource 
allocation preferences among the users. For example, we may consider the weighted sum of rates 

K 

UUR) = J2^^R^ (13) 

i=l 

where Wi > 0, i = 1, . . . , K. When the weights are all equal to unity in (fT3l) . i.e., wi = ■ ■ ■ = wk = 1, 
the utility function is referred to as the sum rate. 

To find the optimal linear precoders, the design variables are the covariance matrices Qi, Q21 ■ ■ ■ , Qk\ 
the precoding matrices Ci are recovered from the covariance matrices according to ©. The optimization 
problem can be formulated as 

maximize U{R) (14) 

over R G R^, G Hf^ (15) 

K- D^i det{IN, + J:^-lH^CkQkClHl') 
subject to Ri < log ^ — - — — " J ' (16) 

" d^^{lN. + Y.ti,k^.H,C,QkClHl^) 

K 

J2ir{EfC,QiCfE,)<P, (17) 

i=l 

where i = 1, . . . , K; j = 1, . . . , B; and R = [Ri . . . Rk]^. Note that within the scope of the optimization 
problem (fT4l)-(fT7l). R, Qi as specified in (fT5l) are simply placeholder variables; they may take on any 
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values in their respective domains as long as all constraints are satisfied. However, the maximization in 
(fT4l)-(fT7]) is not a convex optimization problem; in general, it is difficult to compute the optimal R and 
Qi efficiently. 

IV. Soft Interference Nulling 

A. Precoder Optimization 

In this section, we propose a linear precoding technique called soft interference nulling (SIN), which has 
good performance in the sense that SIN precoding performs better than or equal to any linear precoding 
scheme that aims to completely eliminate interference. The SIN formulation is based on convexifying the 
optimization problem (fT4l)-(fT7l) about some given operating point: Qi G H^^', i = 1, . . . , K . Specifically, 
for the inequality constraints in (fT6l) . we consider the first-order Taylor series expansion of the log- 
determinant function about Qi, i = 1, . . . , K 

K K 

Ri = logdet(^Ir,, + J2HiCkQkClH^)-\ogdet(^I^^+ HiC^QkClH^) (18) 

k=l k=l,k=/=i 
K K 

>\ogdet{lN^ + Y,H,CuQkClH^)- J2 tv{Yr'H,CkQkClH^) 

-logdety,+ tr{Yr^H,CkQkClHf) 

= Ri (20) 
where Yi represents the estimated covariance of the interference plus noise 

K 

Y,^In^+ J2 HiCkQkC^H, eH^^V. (21) 

k=l,kjti 

Note that when Qj's are near Qi's, Ri's are well approximated by Ri's: equality holds when Qi = Qi, i = 
1, . . . , K. Moreover, the inequality in (fT9l ) is due to trX being a global over-estimator of logdet(/„ + X), 
which follows from the concavity of the log-determinant function [|27L 

We consider the precoding matrices that correspond to the solution of following optimization problem 

maximize U (R) (22) 
over ^ e R^, Qi e (23) 

K K 



subject to R^ < logdet(/iv, + J2HiCkQkC^H^) - ^^iYr'H.CkQkCl H,'' 

k=l k=l,k ■ 

K 

-logdety,+ Y tr{Yr'HiCkQkClW^ 



k=l k=l,kj^i 
K 

Ll l^l. lli'^k^k'^f^ ll^^ 

k=l,k^i 

K 

Y,^v{E^CiQiCjE^)<P, (25) 

i=l 

where i = 1, . . . , K; j = 1, . . . , B; the constant matrices Ci, Ej are as given in Q, (fTOl) : and Yi is as 
defined in (|2T]) . Note that the maximization (l22l)-(|25l) above is a convex optimization problem, and its 
solution can be efficiently computed using standard convex optimization numerical techniques, e.g., by the 
interior-point method ||27l , [|28l . The number of active spatial streams for User i, which corresponds to the 
rank of its input covariance matrix Qi, is determined by the optimization framework. In the optimization 
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formulation above, though there is no explicit rank constraint on Qi, we do not expect the effective rank 
of its solution Q* be greater than Ni, i.e., a user cannot receive more spatial streams than it has number 
of receive antennas. 

In the following, we consider an iterative algorithm, as listed below in Algorithm [T] First, we initialize Qi 
to be zero matrices, i = 1, . . . ,K, and solve (|22l) - (|25l) . Li the initial iteration, (fT9l) is a good approximation 
for small interference levels. In the next iteration, each Qi takes on the value of Q* (Algorithm [H 
Line |7]), the solution of the optimization problem (|22l) -(|25]). As we apply the first-order Taylor series 
expansion about the updated Qi in the new iteration, accordingly, the formulation in (fT9l) becomes a good 
approximation about the updated interference levels. This procedure is repeated to iteratively refine the 
estimated covariance Yi of the realized interference plus noise in (|2T]) . Subsequently, we show that the 
iterations are monotonically nondecreasing, and we terminate the algorithm if the utility improvement 
is less than a given tolerance £ > 0. Therefore, when the algorithm terminates, the achievable rates Ri 
(fTSi) are locally well-approximated (to the first order) by the convexified rate expressions Ri in (|20|) 
about the realized interference levels. At the conclusion of the algorithm, the realized interference-plus- 
noise covariance is represented by Yi, there is no stipulation on the interference levels being small. The 
approximation formulation in (fT9l ) is valid for arbitrary numbers of transmit and receive antennas, and 
arbitrary levels of interference in the network. 



Algorithm 1 Iterative Precoder Optimization 

1: Initialize: R <- Okxi, Qi ^ Om^^m^, i = l,...,K 

2: loop 

3: Compute R*,Ql, i = 1,. . .,K, by solving (I22l)-(l25]) 
4: if U{R*) - U{R) < 6 then 
5: break 
6: else 

7: R^R*, Qi^Q*, t = l,...,K 

8: end if 

9: end loop 



We now show that the SIN precoding algorithm performs at least as well as any linear precoding 
scheme that completely eliminates all interference. Recall from Section IIII-B[ a linear precoding scheme 
is interference-free if 

K 

^ ^ HijGjkUk = Ojv^xi, i = l,...,K. (26) 

For example, zero-forcing beamforming achieves the above interference-free condition. 

Proposition 1. The SIN precoding scheme given in Algorithm \T\ performs better than or equal to any 
linear interference-free precoding scheme. 

Proof: Suppose a linear interference-free precoding scheme has precoding matrices: Qi, i = 1, . . . , K. 
Then User i achieves the rate 

Ri = log det (In, + HidQiCjHf) (27) 

where (l27l) follows from substituting the interference-free condition (|26l) in (|7]). Note that in the initial 
iteration in Algorithm[Il Qi = ^MixMi^ and hence Yi = OwixNi, i = I, ■ ■ ■ , K. Substituting (|26|) . (l27l) in 
(|24|) . it is seen that R, Qi are feasible in the SIN optimization problem (I22l)-(l25l). Naturally, the solution 
of the optimization problem is better than or equal to any feasible solution. Furthermore, since Ri is 
a global under-estimator of Ri, the objective function can only improve (or stay the same) after each 
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Fig. 2. Line network model. 
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iteration. As the user rates are bounded, the algorithm converges to a local optimum. Consequently, when 
the algorithm terminates, the SIN precoder performs at least as well as the interference-free precoder. ■ 
An interference-free precoding scheme can be interpreted as one that imposes an infinite penalty on 
interference, whereas SIN relaxes such restriction and allows the possibility of nonzero interference. 
Moreover, SIN precoding is well-defined even when the number of transmit antennas is less than the total 
number of receive antennas. A separate user selection step is not necessary: under SIN precoding, the set 
of active users (i.e., those with nonzero rates) is determined by Algorithm [T] 

B. Clustering Algorithms 

In Section IIII-A[ we assume the coordination cluster Bi for each User i is given. As exhaustive search 
for the best clustering combinations has complexity too high even for small network sizes, in this section, 
we consider two simple clustering techniques. The users choose their coordination clusters based on the 
long-term channel conditions (i.e., the shadowing realizations), but they cannot adapt the coordination 
clusters based on the fast Rayleigh fading realizations. In each algorithm, we assume the coordination 
cluster size \Bi\ > 1 for each User i is given. 

a) Nearest Bases Clustering: In this simple clustering scheme, each User i chooses the nearest (in 
terms of signal strength) \Bi\ bases to comprise its coordination cluster. 

b) Nearest Interferers Clustering: In the nearest interferers clustering scheme, each user attempts to 
reduce interference to its — 1 closest neighbors. First, each User i chooses the nearest base (in terms 
of signal strength) as its home base station. Let B^ denote the home base station of User i. Next, let Xj 
denote the set of neighbors of User i who suffer the most interference from Bi. User i then forms its 
coordination cluster as Bi |J {Bk \ k E Xi}. 

V. Numerical Results 

A. Line Network 

1 ) Network Geometry: We first consider a simple line network model to gain intuition on the behavior 
of different downlink precoding algorithms. Suppose there are B base stations in the network, and they 
are positioned along a line with distance apart. Each base station serves one mobile user, and we 
assume each user is located at a distance dy away from its base station, as illustrated in Fig. [21 The total 
number of users in the network thus is K = B. To minimize the boundary effects, we consider a line 
network with wraparound where the distance dij between User i and Base j is given by 

d,,^\ldl+{dj{t,j))\ t,j = l,...,B (28) 

where 

d{i,j)=dE{0,l,...,B-l} \ d = i-j (mods). (29) 

The radio signal propagation from Base j to User i is modeled as independent Rayleigh fading with 
a distance-based power attenuation of d^^^, where r] is the path loss exponent: i.e., each channel entry is 
independent and identically distributed (i.i.d.) as CJ\f{0, dj^). We assume the line network has geometry: 
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dx = dy = 1, and = 4. In the line network, we consider the case where each base station and each user 
has a single antenna, i.e., M, = Ni = 1, for all i = 1, . . . , K, and j = 1, . . . , B. For such single-input 
single-output (SISO) channels, the channel gain is a scalar, and we will use the notation 

K.^H^^eC, ^ = l,...,K, j = l,...,5. (30) 

2) Dirty Paper Coding: For collaborative cellular networks, a performance upper bound is obtained 
when all base stations cooperate perfectly. Suppose each base station knows the messages of all users, and 
we allow joint encoding at the base stations. Then this cooperative cellular system may be modeled as a 
broadcast channel (BC) with K single-antenna receivers, and an i?-antenna transmitter under per-antenna 
power constraints (PAPC). For a Gaussian MIMO BC, its capacity region [11 J is achieved by the dirty 
paper coding (DPC) scheme [29 1. In DPC, the messages for the users are encoded in a given order, and 
the interference from the previously encoded users is pre-subtracted at the transmitter for the subsequently 
encoded users. In [[221, it is shown that the sum rate of a MIMO BC under PAPC can be computed by 
solving the following convex minimax optimization problem: 

K B 

min max log det Sj/ij/if^ + diag(g) j — log gj (31) 

^ 1=1 i=l 

over g G E^, s G R+ (32) 

K B 

subject to ^ Si < ^ Pj (33) 

i=i j=i 

B B 

i=l j=l 

where 

k = [ha . . . KbY e C^, g = [gi . . . qbY e R?, s^[s,... skY e R^. (35) 

3) Full-Network Zero-Forcing: In the case where all base stations participate in the joint encoding of 
the user messages, a simple linear transmit precoding technique is zero-forcing (ZF) beamforming [9|. 
Unlike DPC, no interference pre- subtraction is performed at the transmitter. Instead, the transmitter sets 
the precoder to be the pseudoinverse of the channel matrix, such that interference is zeroed out at each 
mobile user. Zero-forcing in MIMO BC subject to PAPC is considered in [[T9l - [[2Tll . In particular, the ZF 
sum rate under PAPC can be found by solving the following convex optimization problem: 

K 

maximize E^ log(l + li) (36) 

i=l 

over 7 G R^ (37) 
subject to \W\^-i <[Pi ... PbY (38) 

where 7 = [71 . . .'^kY ^ ^ = [^u]^' the constraint in (l38l) represents component- wise inequality; 
and \W\'^ denotes the component- wise squared magnitude of the entries of W , i.e., = In 
general, ZF is suboptimal. However, at high SNR with multiuser diversity, ZF performs asymptotically 
close to the optimal DPC scheme JH-im. 

4) Achievable Rates and Capacity Bounds: We consider the sum rate of the downlink channel, normal- 
ized by the number of base stations. All the base stations are under the same transmit power constraints: 
Pi = ■ ■ ■ = Pb = P, and we also refer to P as the SNR of the system. In the numerical experiments, 100 
sets of random channel realizations are generated. For each set of channel realizations, the non-cooperative 
and the cooperative normalized sum rates are calculated. Then the rates are averaged over the random 
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Fig. 3. Soft interference nulling (SIN) with different coordination cluster sizes (labeled next to their corresponding curves) and full-network 
zero-forcing (ZF) rates in a line network (B — K = 21, dx = dy — 1), as compared to full-network DPC and non-cooperative full-power 
transmission. 



channel realizations. The convex optimization problems are solved using the software package SDPT3 

moj. 

The line network achievable rates and capacity bounds are shown in Fig. |3] as a function of the SNR 
P, with B = K = 21. In the line network, for simplicity, all users are assumed active and the SIN 
rates are computed with a single iteration of Algorithm [TJ The non-cooperative baseline refers to the 
case where each base station transmits at full power. Without base station cooperation, the average user 
rate saturates at approximately 1.6 bps/Hz. On the other hand, interference can be overcome by allowing 
the base stations to cooperate, as shown by the DPC rates when the cooperative system is modeled as 
a BC. ZF beamforming is able to achieve increasing system throughput as the SNR improves. There is 
a gap between between the ZF rate and the DPC rate, but they both exhibit similar scaling trends as P 
increases. In particular, the ZF beamforming technique allows the network to overcome its interference- 
limited performance bottleneck, thus demonstrating the value of cooperative cellular networks. However, 
ZF beamforming requires all base stations in the network to cooperate. 

5) Soft Interference Nulling Precoding: For the SIN precoding scheme, different coordination cluster 
sizes are considered, and the cluster sizes are labeled next to their corresponding curves on the plot. For 
the line network topology, the two clustering algorithms discussed in Section IIV-BI are equivalent: i.e., 
each user chooses the closest \Bi\ base stations to participate in its coordination cluster. For example, when 
the coordination cluster size is 3, User 1 is served by Bases {21, 1, 2}, User 2 by Bases {1, 2, 3}, User 3 
by Bases {2,3,4}, and so on. A coordination cluster size of 21 represents full network coordination. 

Note that the full-network ZF beamforming scheme satisfies the interference-free condition (|26l) ; there- 
fore, under full network coordination, by Proposition [H SIN precoding outperforms ZF. In particular, at 
low SNRs, ZF suffers from noise amplification while SIN does not. Moreover, at moderate SNRs, SIN with 
limited coordination cluster sizes is able to outperform ZF with full network coordination. For example, 
at the SNR P = 18 dB, SIN with a coordination cluster size of 7 outperforms ZF with full network 
coordination. As the SNR increases, however, it is observed that SIN under partial network coordination 
once again becomes interference-limited. 
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Fig. 4. Hexagonal three-sectored cellular wireless network. There are 19 cells in the network, with wraparound at the edges. Each cell 
has three sectors. Each sector corresponds to a base station, and each base station serves one mobile. Each arrow represents the boresight 
direction of a base station's antenna beam. 



B. Cellular Network 

1) Network Geometry: We next consider a cellular network, as shown in Fig. HI which consists of two 
rings of hexagonal cells, with wraparound at the boundary. Each cell has three sectors; thus there are 19 
cells, or 57 sectors, in the network. Each sector has one transmit antenna and serves one user, and each 
user has a single receive antenna. Users are randomly populated in the network. The distance between any 
two closest cell centers is 0.5 km. Average channel SNR is determined by propagation path-loss (with a 
path-loss exponent of 3.76) and shadow fading (with 8-dB standard deviation, 50-m correlation distance, 
50% correlation across sites). The transmit antenna at each sector has a parabolic beam pattern. A user 
is associated with the sector to which it has the highest average SNR (up to the maximum of one user 
per sector). The users are indexed such that Base i wishes to transmit to User i. In composition with the 
path-loss and shadowing, each channel also experiences i.i.d. fast Rayleigh fading. Each sector is under 
a transmit power constraint that corresponds to a cell-edge SNR of 20 dB. 

2) Myopic Zero-Forcing: Zero-forcing requires all base stations to cooperate in order to ensure an 
interference-free signal at each of the users. However, in a wide-area cellular network, it is unrealistic 
to demand full-network coordination. When only clusters of base stations cooperate, interference cannot 
be completely eliminated. For comparison, in this section, a simple myopic scheme is considered. First, 
in the case when a base station belongs to multiple clusters, it divides its transmit power equally among 
the clusters. Then, within each cluster, the ZF precoder is used. Finally, the inter-cluster interference is 
treated as noise when the achievable rates are computed. Myopic ZF requires less backhaul bandwidth in 
exchanging CSI; however, [25| shows that the backhaul overhead is dominated by data sharing, for which 
the same bandwidth requirements apply for myopic ZF and SIN. 

3) Performance Comparison: For the numerical experiments, 10 instances of shadow fading realizations 
are generated. For each shadow fading realization, 10 Raleigh fading instances are generated (i.e., a 
total of 100 sets of channel realizations). In the cellular network, we consider a simple user selection 
mechanism. In calculating the myopic ZF rates, the users that achieve the lowest 10% non-cooperative 
rates are allowed to be in outage. To compute the SIN precoding matrices, iterations of Algorithm [T] 
are carried out until a tolerance oi e = 0.01 is achieved (i.e., the set of active users is determined by 
the algorithm). The average rates in cellular network under different clustering algorithms for various 
coordination cluster sizes are shown in Fig. [51 In the figure, clustering (a), (b) refer to nearest bases 
clustering, and nearest interferers clustering, respectively, as described in Section IIV-BI It is observed that 
nearest interferers clustering outperforms nearest bases clustering, which suggests cooperation is more 
useful in mitigating interference than boosting the users' signal strength. Unfortunately, the myopic ZF 
scheme does not satisfy the interference-free condition in (|26l ), and Proposition [T] does not apply in the 
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Fig. 5. Average rates under different transmission strategies for various coordination cluster sizes. Clustering (a) refers to nearest bases 
clustering, and (b) refers to nearest interferers clustering. 



comparison of its peri'ormance with SIN precoding. However, the simulation results show that for a given 
coordination cluster size, SIN precoding achieves considerably higher throughput compared with myopic 
ZF. Furthermore, SIN can more effectively take advantage of larger coordination clusters: the SIN rates 
improve when the coordination clusters become larger; whereas the improvement in the myopic-ZF rates is 
only marginal. Intuitively, SIN precoding is able to achieve better performance over myopic ZF because i) 
SIN optimizes the transmit power allocation among the coordination clusters, and ii) SIN aims to achieve 
the optimal balance between beamforming and interference nulling in using the multiple antennas in each 
cluster. 



VI. Conclusions 

In this paper, we consider maximizing a concave utility function of the user rates in a cooperative 
MIMO cellular downlink network. Without base station cooperation, a cellular network is interference- 
limited. Cooperation among base stations allows the joint encoding of user signals, which can overcome 
the interference limitation. However, the capacity-achieving dirty paper coding (DPC) scheme has high 
complexity. When all base stations in the network cooperate, zero-forcing (ZF) beamforming offers good 
performance relative to DPC. We investigate the case where each user may specify only a subset of the 
base stations to form a coordination cluster, and different coordination clusters may overlap in the network. 
The base stations share CSI, but a user's message is only routed to base stations in the coordination cluster 
associated with the user. 

We focus on linear precoding techniques for low-complexity implementation. In general, the optimal 
linear precoding is a nonconvex problem which is difficult to compute efficiently. We propose an approx- 
imation formulation called soft interference nulling (SIN), which performs as well as or better than ZF 
beamforming. The SIN precoding matrices are computed by solving a sequence of convex optimization 
problems. In a line network, it is shown that SIN precoding under partial network coordination can 
outperform ZF under full network coordination at moderate SNRs. In a hexagonal three- sectored cellular 
network, SIN precoding achieves considerably higher throughput compared to myopic ZF with the same 
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coordination cluster size. Moreover, the SIN rates improve with increasing coordination cluster sizes, 
while the myopic-ZF rates do so only marginally. 
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